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The linear spin-1/2 Heisenberg antiferromagnet with exchanges Ji, J2 between first and sec- 
ond neighbors has a bond-order wave (BOW) phase that starts at the fluid-dimer transition at 
J2/J1 = 0.2411 and is particularly simple at J2/J1 = 1/2. The BOW phase has a doubly degen- 
erate singlet ground state, broken inversion symmetry and a finite energy gap Em to the lowest 
triplet state. The interval 0.4 < J2/J1 < 1.0 has large Em and small finite size corrections. Exact 
solutions are presented up to A*' = 28 spins with either periodic or open boundary conditions and for 
thermodynamics up to A'^ = 18. The elementary excitations of the BOW phase with large Em are 
topological spin-1/2 solitons that separate BOWs with opposite phase in a regular array of spins. 
The molar spin susceptibility xm{T) is exponentially small for T <C Em and increases nearly linearly 
with T to a broad maximum. Ji, J2 spin chains approximate the magnetic properties of the BOW 
phase of Hubbard-type models and provide a starting point for modeling alkali-TCNQ salts. 
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I. INTRODUCTION 



The extended Hubbard model (EHM) has competing 
on-site repulsion U > 0, intersite interaction > and 
electron transfer t between neighbors in one dimension 
(ID) with evenly spaced sites. The half-filled case with 
one electron per site has several phases: a charge den- 
sity wave (CDW) at V > U/2 with broken electron-hole 
symmetry and occupation numbers rt > 1 on one sublat- 
tice, 71 < 1 on the other; a spin fluid phase at = as 
known from Hubbard models; and as proposed by Naka- 
mura, 1] a bond-order wave (BOW) phase with broken 
inversion symmetry between the CDW and the spin fluid 
phases when t/U is sufficiently large for a continuous 
CDW transition. The BOW phase has long-range order 
and a flnite magnetic gap Em to the lowest triplet excited 
state. Multiple theoretical approaches, primarily at 
U < 2t and U — At, have confirmed a narrow BOW phase 
in the EHM. Other spin-independent potentials also sup- 
port a BOW phase when the CDW transition is continu- 
ous . The narrow BOW phase of Hubbard-type models 
presents major computational difficulties. 

In this paper, we consider the BOW phase of a fa- 
miliar spin-1/2 chain with frustrated antiferromagnetic 
(AF) exchange Q- The BOW phase becomes numeri- 
cally accessible and can be demonstrated in finite sys- 
tems. Although charge fluctuations are strictly excluded 
in spin chains, the BOW phase again illustrates broken 
inversion symmetry at sites, long-range order and finite 
Em- The spin chain has AF exchange between first and 



second neighbors. 
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We consider the interval < x < 1 and set the total 
exchange J = 1 as the unit of energy. The x — Q 
limit is a linear Heisenberg antiferromagnet (HAF). 
Second-neighbor exchange J2 = a; J for a; > opposes 
short-range antiferromagnet (AF) order and eventually 
induces a fluid-dimer phase transition that has been the 
focus of recent studies The X ^ \ limit gives 

two HAFs on the even and odd sublattice, respectively. 
White and Affleck [l3| considered J2/J1 > 1 using field 
theory and the density matrix renormalization group 
(DMRG). 

Okamoto and Nomura Q located the transition at 
Xc = — xi) = 0.2411, or Xi = 0.1943 in our no- 

tation, where a magnetic gap Em opens. The "dimer" 
phase refers to the earlier observation of Majumdar and 
Ghosh (MG) that, for an even number N of spins and 
periodic boundary conditions (PBC), the exact ground 
state (gs) at xmg = 1/3 (J2 = ^i/2) has singlet-paired 
spins on adjacent sites, just as in the Kekule diagrams 
shown in Fig. [1] 



\K1) - (l,2)(3,4)...(7V-l,iV) 
\K2) = (2,3)(4,5)...(7V,1) 
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Kekule diagrams are particularly simple BOWs. They 
illustrate broken inversion symmetry, long-range order, 
double degeneracy and finite excitation energies [l3| at 
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FIG. 1: Lines indicate singlet-paired spin on adjecent sites. 



the MG point. Since these are the defining features of a 
BOW phase, that is what we will call the "dimer" phase. 
The notation H{x) is convenient for the x = 1 limit 
of HAFs on the even and odd sublattices. Now x < 1 
describes interchain exchange that is frustrated because 
each spin is coupled to two neighbors of the other 
sublattice. With constant total exchange, the gs en- 
ergy is highest when the spin chain is the most frustrated. 

For open boundary conditions (OBC), \K1) is the ex- 
act nondegenerate gs of (1/3). The chemical analogy is 
now to partial double and single bonds in linear polyenes 
or in polyacetylene. The BOW associated with \K1) is 
well understood in dimerized arrays whose elementary 
excitations are the topological solitons of the uncorre- 
lated Su-Schrieffer-Heeger (SSH) model Similar 
conclusions hold in correlated models of conjugated 
polymers 17[ or ion-radical stacks [18]. Spin solitons 
in BOW systems at finite temperature separate \K1) 



and \K2) regions in infinite regular chains. The Peierls 
instability of Hubbard or spin chains is a separate topic 
that requires electron-phonon or spin-phonon coupling. 
Our discussion of H{x) is limited to regular arrays with 
PBC or OBC. 

To introduce the principal features of H{x), we show in 
Fig. [2]the gs energy per site, eo(a;)jfor intermediate TV « 
20 and PBC. Bonner and Fisher [l^ found that eo(0) of 
the HAF converges as « A^~^ to the exact value, — ln2 + 
1/4, due to Hulthen [2^ and denoted by arrows at a; 
and 1. Convergence at x = 1 is for two HAFs of i V/2 
sites, from below when A^/2 is even and from above |21| 
when N/2 is odd. The shape of eo(a;) indicates different 
frustration at small x for exchange J2 in one HAF and at 
large x for exchange Ji between two HAFs. Frustration is 
greatest at the eo (x) maximum. The Hellmann-Feynman 
theorem gives 
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FIG. 2: Ground state energy per site, £0(2:), of the spin chain 
H{x) in Eq. [l]for A'^ sites with periodic boundary conditions 
(PBC). The BOW phase starts at the fluid-dimer transition at 
xi; X2 is the excited state crossover discussed in the text. The 
energy of \K1) or \K2) in Eq. [2] is exact at xmg = 1/3 and 
approximate for x < 1/2; the energy of a product of Kekule 
diagrams \KA)\KB) is approximate for x > 1/2. Arrows at 
a; = and 1 mark eo of the infinite chain. 



at a; = than at a; = 1. First and second neighbor spin 
correlation functions are known [22j] exactly at x — 
and they add in Eq. [3] The slope at a: = 1 is just the 
first-neighbor correlation function. 

The solid lines in Fig. [2] are the energy of \K1) or 
\K2) for X < 1/2, exact at xmg = l/3j and of a product 
of Kekule diagrams of two HAFs for x > 1/2. The BOW 
phase that we characterize below starts at xi. We found 
X2 ~ 2/3 using the Okamoto-Nomura [1] treatment of 
xi. The gap is exponentially small but finite for 
X > X2, and the BOW phase extends to x = 1 according 
to White and Affleck Q. 

The order parameter B{x) is the gs amplitude of the 
BOW, 
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The two gs have ±B{x). It follows immediately that 
5(1/3) = 3/8 for \K1) or \K2) in Eq. H As shown 
below, large B{x) and E„i{x) between x = 1/3 and 
X « 1/2 make possible our detailed finite-N study of the 
BOW phase. 



The bond orders or spin correlation functions are equal 
at deo /dx = when there is equal choice for pairing with 
a first or second neighbor. The slope de^/dx is steeper 



The complete basis of H{x) has dimension 2 , since 
each spin-1/2 has two orientations, and the total spin 
< 5 < N/2 is conserved. Reflection a through sites 
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corresponds to inversion symmetry Ci at sites in the infi- 
nite chain. Valence bond (VB) methods [23l - [26| are well 
suited for finite models that conserve S. A few states 
with any S and a can be found exactly up to TV « 30. 
The full spectrum is needed for thermodynamics and has 
been obtained [l^] to iV = 16, which we increase to 
TV = 18. DMRG extends ^ thermodynamics to N = 
64. The spin chain H{x) benefits from the smaller basis 
compared to 4^ in Hubbard models with charge degrees 
of freedom. An even greater advantage may be the exact 
gs at the MG point for finite A'^. In contrast to the nu- 
merically difficult BOW phase of Hubbard-type models, 
the BOW phase of H{x) is accessible to direct finite-A 
modeling between a; w 1/3 and x w 1/2. 

We characterize the gs properties of the BOW phase 
and its elementary excitations in Section II, including the 
magnetic gap Em{x), the order parameter B{x), excited 
states at the MG point, and the bond-order domain walls 
of spin solitons. The temperature dependence of the mo- 
lar spin susceptibility xm{T) and specific heat C(T) are 
found in Second III. Following an activated regime that 
depends of xm(2^) increases almost linearly with T 
in the BOW phase, quite differently from an HAF or an 
EHM with t<^{U-V). The Discussion relates H{x) to 
the EHM and to vr-radical salts with xm (T) nearly linear 
in T. 



II. GROUND AND LOW-ENERGY STATES 

We use valence bond (VB) methods [2^, [2^ to solve 
H{x) exactly for finite A" and either periodic or open 
boundary conditions in exact subspaces with fixed total 
S and reflection cr = ±1 at sites. The gs is a singlet, 
5" = 0, with either a = 1 or -1 depending on A^ and x. 
The linear combinations |A'l) ± |Ar2) in Eq. [2]transform 
as a = ±1, even or odd under inversion in the infinite 
chain. The gs for other a: is a symmetry adapted linear 
combination of singlet VB diagrams |fc). We define E„{x) 
as the excitation energy to the lowest singlet with oppo- 
site a symmetry. Figure [3] compares E„ {x) for A = 24 
and PBC to the gap Em{x) to the lowest triplet. Finite- 
size effects are large at a; = for an HAF of 24 sites 
and about twice as large at a; = 1 for two HAFs of 12 
sites, as expected when excitations energies go as « 1/A^. 
We obtained similar graphs of E„{x) and Em(x) up to 
N = 28. There is no difference at small x aside from 
1/A^ effects. The excitations are qualitatively different 
at a; w 1, however, when A^/2 is even or odd. Even A/2 
is required for proper comparison at a; = and 1. When 
A^/2 is odd, the a; = 1 limit corresponds to two HAFs 
with an odd number of spins and a doublet gs. Since the 
radical also has two- fold orbital degeneracy |2lj, there 
are several gapless excitations in the x = 1 limit of no 
interchain exchange. 

Okamoto and Nomura [sj identified the quantum 
transition at xi by finding Ea-{xi) = Em{xi) from A^ 




0.2 0.4 0.6 0.8 

X 



FIG. 3: Finite-size effects on excitation energies of H{x) for 
N = 2A: with PBC. Em is the excitation to the lowest triplet 
and Ecj to the lowest singlet with opposite inversion symme- 
try. The crossovers xi{2A) and X2{2A) are listed in Table U 



= 10 to 24. They argued that an excitation crossover 
at finite A^ is more accurate than extrapolation to find 
Em{xi) — 0. The slow variation of a;i(A^) and extrapola- 
tions made possible their accurate determination of Xi. 
Our results for xi{N) in Table 1 to A^ = 24 agree with 
ref. 8. Previous work 0-03 focused on the fluid-dimer 
transition at xi, while we are interested in the BOW 
phase with x > xi. The same method yields Em = Ea 
at X2{N) in Table U for even and odd 7V/2. Finite-size 
effects are stronger because the chains are effectively 
half as long. A joint 1/A^ extrapolation of the two 
sequences returns X2 = 0.67±0.01 (J2/J1 = 2.03±0.03). 
The X2 crossover does not signify the termination of 
BOW phase, however, which extends [lii] to x = 1. We 
improved the accuracy of the DMRG algorithm [2^ to 
look aX X > X2 and find small but finite Em{x) and B{x) 
up to a; = 0.8 (J2/J1 = 4.0), beyond which even more 
accurate DMRG is required. We do not understand 
the different implication of xi and X2 crossovers, but 
note that H{x) also has a spiral phase (s^l starting at 
xmg = 1/3 whose order parameter is twist angle. In 
the present work we focus on the BOW phase with large 
B{x) and Em{x). 

The BOW phase of the extended system has degener- 
ate gs in the a — ±1 sectors and hence Ea{x) — 0. Finite- 
size effects are extraordinarily small for 0.3 < a; < 0.5 
where Eq is close to —N/A. The difference in total en- 
ergy, Eo{x, 1) — Eq{x, — 1), for cr = ±1 is shown in Fig. 0] 
up to A^ 30. The absolute gs for x < 1/3 is £;o(-l) for 
N ^Ap and Eq{1) for A^ = 4p 2. The Eo{x, ±1) curves 
are degenerate at x = 1/3 without crossing. They cross 
for a; > 1/3 and the number of crossings depends on A^. 
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TABLE I: Crossing points xi{N) and X2{N) where Em{x) 
Ea{x) for iV sites and PBC. 



N 




X2 {N/2 odd ) 


X2 {N/2 even) 


18 


0.1949 


0.5669 




20 


0.1947 




0.6262 


22 


0.1947 


0.5784 




24 


0.1946 




0.6368 


26 


0.1946 


0.5885 




28 


0.1944 




0.6421 


oo 


0.1943 " 


0.68 


0.66 
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FIG. 4: Finite-size effects on the energy the lowest singlets of 
H{x) with opposite inversion symmetry in the BOW phase. 
The energies are equal at a; = 1/3 without crossing. Crossings 
at a; > 1/3 are governed by the ground-state symmetry at 
a; = and 1 as discussed in the text. 



When N — Ap and p is even, the gs at a; = 1 also trans- 
forms as a = —1; there is an even number of crossings 
for iV = 24 or 16. The E„ bump in Fig. [3] at x = 0.5 is 
due to two crossings. For N — Ap and odd p, the x = 1 
gs has (7=1 symmetry that requires an odd number of 
crossings for = 28 or 20. The gs for = 4p + 2 has 
a = 1 symmetry for x < 1/3. There is an even number 
of crossings up to a; = 0.60 for = 26 or 18 (even p) and 
an odd number for N = 30 or 22 (odd p). Degeneracy 
without crossing at the MG point and subsequent sym- 
metry crossovers for finite N are the principal reasons for 
remarkably small E„{x) in this interval. 

Since Ea{N) > is due to finite-size effects in the 
BOW phase, extrapolation of Em{N) — Ea{N) yields 
Em{x). As seen in Fig. [SJ E„i{N) — E„{N) converges 
well to Em{x) on the xi side and less well on the X2 side 
where even and odd A^/2 appear for 3:2 (iV) in Table HI 
The largest magnetic gap is 0.29J at x = 0.40, close to 
the eo(a;) peak in Fig. [2] and clearly beyond xmg = 1/3. 
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FIG. 5: Excitation energy Em{x) to the lowest triplet of 
H{x) for finite A'' and PBC in the BOW phase. The infinite 
chain has Ec{x) = and a doubly degenerate ground state. 



The Em maximum and position agree well with DMRG 
in Fig. 5 of ref. 14. Since Fig. 2] shows Ea-{x) to be very 
small between x = 0.3 and 0.5, large Em in this inter- 
val is consistent with small finite-size corrections. More 
accurate DMRG is needed [H] for Em{x) at a; > 0.5. 

To obtain the BOW amplitude B{x), we break inver- 
sion symmetry according to 



H{x,6) = H{x) 



(5) 



At the MG point, \K1) is the gs for (5 = 0-f and 
\K2) for 5 ~ 0—. The gs energy per site, eo{x,S), 
gives B{x) — —{deo/dS)o. The inset of Fig. [HI shows 
-(e(l/3,(5) - e(l/3,0))/(5 for A^ = 20 as a function of 
S. The intercept is B{l/3) = 3/8 while the slope is 
Xd/2, the harmonic electronic force constant per site for 
dimerization that will be needed in a later study of lat- 
tice vibrations. Fig. [6] shows B{x) in the BOW phase. 
5(1/3) = 3/8 follows directly from — (s„.s„+i) = 3/4 or 
for paired and unpaired neighbors, respectively. The N 
dependence of B{x) is negligible near the MG point up to 
X ~ 0.45, but it becomes significant around a; « 0.5 where 
the location of the gs crossings in Fig. |4] depend on A'^. 
Finite-size effects also appear near Xi and X2 where B{x) 
becomes small but does not vanish. Since ^/'o(a;',±) with 
CT = ±1 are degenerate at crossings x', the linear combi- 
nations ("00 {x' , +) ± ipo{x' , — )) / a/2 are broken-symmetry 
states whose expectation value in Eq. |4] leads to 



B{x') = |(^o(a:',+)|^(-l)"s„.s„+i|V'o(x',-))|/A^ (6) 
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FIG. 6: Amplitude B{x) in Eq. |4]of the BOW of H{x) for 
finite N and PBC in the BOW phase. B(l/3) = 3/8 is exact 
for either \K1) or \K2) in Eq. [2] Finite-size effects around 
X ~ 1/2 are due to reversals of the ground states inversion 
symmetry. The inset shows — [eo(l/3, 5) — eo{l/3)]/5 vs S at 
A'^ = 20 for the symmetry-breaking perturbation in Eq. [S] 



The matrix element agrees quantitatively with B{x') — 
— {deo{x' ,6)/d5)Q, as it must. The two determinations of 
B{x) are the same within our numerical accuracy when 
E„ < O.OIJ. 

To our surprise, -8(0.35) is slightly but distinctly 
larger than 3/8, the amplitude at the MG point for \K1) 
or \K2). A Kekule diagram has perfect AF correlation 
with one neighbor, which seems to be the limiting case 
of a BOW. While the AF correlation or bond order 
decreases slightly in systems with B{x) > 3/8, there is 
now small F correlation or negative bond order with the 
other neighbor. Direct solution of systems with OBC 
and B{x) > 3/8 yields large positive and small negative 
bond orders that alternate along the chain. We recall 
that the second-neighbor bond orders, — (s„.s„+2), are 
negative for x < 1/3, vanish at a; = 1/3, and are positive 
for X > 1/3. The BOW phase for a; > 1/3 has AF 
correlations for second neighbors and alternating strong 
AF and weak F correlations for first neighbors. We also 
note that B{x) and E„i{x) are not simply proportional 
to each other. The B{x) maximum in Fig. [5] is at 
decisively lower x than the Em{x) maximum in Fig. [5] 
White and Affleck 14] were also surprised that d = 2B 
could exceed 3/4 at J2/J1 > 1/2 and interpreted the 
result as ferromagnetic correlation; Fig. 8 of ref. 14 
agrees quantitatively with Fig. |6]for 0.3 < a; < 0.5. 

We consider next the excited states of H{x) and 
present results at xmg = 1/3 that are representative 
for the interval 0.3 < x < 0.5 in which B{x) and Em{x) 
are large. The gs energy per site at the MG point is 



TABLE II: Excitation energies E{Sr) of H{l/3), in units of 
J, for N sites and PBC. 

Spin and State, Sr N=28 N=26 N=24 



Triplet, Tl 


0.1691 


0.1705 


0.1727 


Singlet, S3 


0.1757 0.1793 


0.1839 


T2, T3 


0.1833 


0.1873 


0.1921 


S4, S5 


0.1898 


0.1953 


0.2022 


T4 


0.2089 


0.2155 


0.2242 


Quintet, Ql 


0.4200 


0.4324 


0.4499 



eo(l/3) = -1/4 for either PBC or OBC. Table H] list ex- 
citations with increasin g en er g y for N — 24, 26 and 28. 
Sparse matrix methods [25|, [2a] are used for a few states 
in each symmetry subspace. It becomes progressively 
more difficult numerically to go beyond 3 or 4 states for 
large N. The notation Sr indicates total spin S and state 
index, r = 1,2,3.... States are doubly degenerate with 
wavevector ±fc except for k — 0(cr = +1) and 7r(fT = — 1). 
The lowest triplet at Em and singlet at E3 decrease slowly 
with N and are known rigorously to be finite in the in- 
finite chain [T^. Finite-size effects are more pronounced 
with increasing r. The gap E3 — Em decreases with N. 
We expect it to vanish in the extended system whose el- 
ementary excitation are spin solitons, each with s = 1/2, 
with paired or parallel spins. There are additional sin- 
glets and triplets below the lowest quintet at Eq{N). 

Since all sites of H{x) are equivalent for PBC, it is 
difficult to discern solitons even with exact eigenstates 
in hand. Fortunately, the gs energy per site for OBC 
and even N is again cq ~ —1/4, and the gs \K1) in 
Fig. [1] has alternating bond orders of 3/4 and along 
the chain. We consider -ff(l/3) with OBC and odd N, 
either = 4p -t- 1 or ./V = 4p — 1. The gs is a doublet, 
S — = 1/2, with spin density pn = 2(5^) at site n. 
With central site at n = 0, the terminal sites are ±2p 
when N ^ Ap + 1 and ±2{p - 1) when N = Ap - 1. 
Linear polyenes or VB diagrams rationalize two distinct 
series when N is finite. The pentyl radical {N — 5) 
has po > and partial single bonds at the center, while 
the allyl radical {N = 3) has po < and partial double 
bonds at the center. Soliton spin densities of -ff(l/3) for 
N = 25 and 23 are shown in Fig. [71 Sites with p < 
indicate electronic correlation |3l| and correspond to 
nodes in uncorrelated Hiickel or tight-binding theory. 

The gs bond orders of -ff(l/3) for odd N are close 
to 3/4 at the end and reverse in between. The {N — 
l)/2 bond orders are symmetric about the center, n — 
0. Fig. [8] displays — (s„.s„+i) for different values of N. 
Bond orders oscillate with increasing n and grow from 
the center. As expected, the central bond order is slightly 
larger for the 4p — 1 series than for the ip + l series. Both 
spin densities and bond orders are typical of spin-1/2 
solitons that for H{l/3) connect \K1) and \K2) regions. 
In a BOW phase, broken inversion symmetry and solitons 
are found in regular arrays. Of course, the soliton width 
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2^ depends on models and parameters; 2^ increases with 
decreasing dimerization in the SSH model [15| and it also 
depends on correlations. The results in Figs. [7] and [8] 
suggest that spin solitons at the MG point have 2^ « 15. 
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FIG. 8: Ground-state bond orders, — {sV,.s„+i) of ff(l/3), 
for odd N and open boundary conditions (OBG). The bond 
orders are symmetric about the center, n — 0, and increase 
to almost 3/4 at the ends. 



FIG. 7: Spin densities p„ = 2(5^) in the doublet ground state 
of H{l/3) with OBC and N = 23 and 25. The central site 
at n = has large positive po for A'' = 25 and small negative 
po for N = 23, as discussed in the text. The spin soliton is 
delocalized over the central part in either case. 



The energy 2Ew of two domain walls is found by com- 
paring the gs energy of even and odd systems with OBC 
and equal length, 

2Ew{x, N) = 2Eo{x, N) - Eo{x, N 1) - Eo{x, N + 1)(7) 

At the MG point, we find 2Ew= 0.1701, 0.1684 and 
0.1669 for N =23, 25 and 27, respectively, slightly less 
than the S„ values in TableHIfor even N. 2Ewil/5, N) 
has weaker N dependence than £'m(l/3, TV), and a joint 
extrapolation returns 2Eyy = Em — 0.151 for the infinite 
chain. Finite-size effects are larger for i?3(l/3,iV) in 
Table In] and even larger for Eq{1/3,N). 

Direct solution up to ~ 30 indicates that the ele- 
mentary excitations of H{l/3) are spin-1/2 solitions with 
2Eyj = Em — Es = Eq/2 in the infinite chain. Finite N 
results suffice for 0.3 < a; < 0.5 when B{x) and Em{x) are 
large. Longer chains can be studied using DMRG meth- 
ods that will be needed for the BOW phase of Hubbard- 
type models. Since ID systems at T > cannot have 
long-range order, topological solitons are generic features 
of systems with a BOW phase. The present discussion 
is limited to a rigid lattice with purely electronic domain 
walls, but solitons are also expected in deformable lat- 
tices with linear electron- or spin-phonon coupling. 



III. MAGNETIC SUSCEPTIBILITY AND 
SPECIFIC HEAT 

Static magnetic susceptibility provides by far the most 
direct comparison with experiment, as amply illustrated 
[3^ [33} by Heisenberg and other spin chains and by spin- 
Peierls systems. The molar spin susceptibility, XAiiT), is 
an absolute comparison for organic radicals with small 
spin-orbit coupling and g close to 2.00236, the free- 
electron value. Since H(x) conserves S, the energy level 
Esr splits into 2S + 1 Zeeman levels in an applied field. 
The full energy spectrum of H{x) in zero field is required 
to construct the partition function 

N/2 

Qn = Y1 + l)eM-Esr/kBT) (8) 

S=0 r 

where /cs is the Boltzmann constant and Esr is excitation 
energy from the singlet gs to the r*'' energy level with 
spin S. The molar spin susceptibility of an iV-site chain 
is [13 

22 

i2S+l)eM-Esr/kBT) (9) 

where Na is Avogadro's number and /is is the Bohr mag- 
neton. Finite-size effects become severe when ksT is 
small compared to Em{N). 

Fig. inishows xm for N = 16 as a function of ksT/ J for 
several values of x. The curves converge for T > J/ks 
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FIG. 9: Molar spin susceptibility, xm{T,N) in Eq. (10), vs 
ksT/J for spin chains with iV = 16 and PBC for x = Q, 0.15, 
0.25 and 0.40 in Eq. (1). The fluid phase up to xi = 0.1943 
has finite x m (0) . The curves become independent of x at high 
T. 



FIG. 10: Molar susceptibility, Xm{T, N) in Eq. El vs ksT/Ji, 
the nearest-neighbor exchange, for spin chains with A'^ = 18 
and PBC. The dashed lines for x = 0.25 and 0.40 are Eq. [11] 
with parameter shown in Fig. 1111 



because the total number of spins is the same and so 
is the Weiss constant J/2fcs- The number of spins and 
Weiss constant are the and terms, respectively, 
at high T. The XMiT) maxima are well converged at 
N = 16, as can be shown by solving iV ~ 18 or 14. 
The situation is different as T — >■ 0, where x < xi leads 
to finite xm(0) in the fluid phase Q while x > xi has 
Xm(0) = due to finite Em- The x = and 0.15 curves 
for finite N are dominated by finite-size effects at low 
enough T. The x = 0.40 curve is almost quantitative 
since Em{OA) exceeds zero-field effects at = 16. The 
X — 0.25 curve is intermediate since Em is finite but N 
dependent. 

For reasons given in the Discussion, we are interested 
in variable J2 = xJ at constant Ji = J(l — x). The 
Xm(T)Ji maxima in Fig. [TOl depend weakly on J2 up 
to a: = 0.40; the curves now cross because the Weiss 
constant varies with x. The HAF {x = 0) maximum 
broadens and shifts to lower T with increasing Jij J\ = 
xj (1 — a;). Size convergence at low T is shown in Fig. [TT] 
for X = 0.25(J2/Ji = 1/3) and x = 0.40 (J2/J1 = 2/3). 
Large E^ at a; = 0.40 gives convergence at A^=16 and 
18. Small Em at x = 0.25 limits convergence to the broad 
Xm(T) maximum. In either case, finite Em ensures that 
Xm (0) = and gives a substantial range in which xm (T) 
is almost linear in T. The slope of the linear regime 
depends weakly on N . 

There are many realizations of dimerized HAFs with 
a; = in Eg. [T]and alternating J(l ± 5) along the spin 
chain [s^ls^. Dimerized chains with 5 > 0.3 are well ap- 
proximated as A T/ 2 singlet-triplet (ST) pairs with spin- 
wave dispersion (s^ . Issj . An ST approximation also pro- 
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FIG. 11: Finite-size effects on xm{T,N) in the BOW phase 
of H{x) at X = 0.25 and 0.40. Large i5m(0.40) gives small 
changes for N = 16 and 18, while small iJm(0.25) leads to 
stronger A'^ dependence. The dashed lines for x = 0.25 and 
0.40 are Eq. [TT]with the indicated parameters Zi and A. 



vides insight into the BOW phase. A normalized density 
n{E) of two-level systems with ST gap E leads to 



Xst{T) 



n{E)dE 



3 + c^p{E/kBT) 



(10) 
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The integral can be evaluated for any piecewise constant 
n{E). We consider n{E) — aside from two intervals; 
n{E) = A/ {El - E^) for E^ < E < Ei and n{E) = (1 - 
A)/{E2-Ei) for Ei<E<E2. The molar susceptibility 
is 



maximum shifts smoothly to higher fc^T/Ji in the fluid 
phase with a; < xi. Likewise, there is a single C(T) peak 
when E„i becomes small for x > 0.5. Finite E„i in the 
BOW phase shifts the singlet and triplet part of C(T) to 
low energy and separates them from S >2 contributions. 



Xst{T, N).h 



A 



3(zi 
1 - 



-In 



A 



1 + 3exp{-ZjnJikBT) \ 
1 + 3exp(— zi Ji/ssT) / 



3(^1 - Z2) 



In 



1 + 3exp(— zi JifcsT) 
1 + 3exp(— Z2 Jifc^T) 



(11) 



with z„i = E„iJi, zi = El J I and Z2 — £'2>/i- As 
expected, the ST gap gives an exponential xst sX 
sufficiently low T, thereby fixing Zm- The width of 
n{E) is controlled by Z2 and is fixed by the Xm{T) 
maximum. The shape of xst{T) can be varied by A 
and zi, or by additional parameters when n[E) has 
more than two intervals. The dashed lines in Figs. [10] 
and [11] are xst [T) with the parameters in Fig [TT] The 
spin susceptibility in the BOW phase is reasonably well 
modeled with a distribution n{E) of ST gaps that is 
constant in two intervals. 

The molar specific heat of H{x) for N sites is 



C(T) _ 1 / J 
NAkB ^N\k^ 



{{EiTf) - {E{T)f) (12) 



The thermal averages require the energy spectrum Esr 
and degeneracy. The results below are for PBC. Since S 
is conserved, separate contributions to {E{T)) can read- 
ily be identified. But the entropy is not additive in S 
and there is no unique partitioning of C(T)dT — TdS 
into contributions in S. One choice is the temperature 
derivative of the S component of {E{T)). Another choice 
is based on fluctuations, 



C{T) 

NAkB 



NQ 



N 



knT 



)'^(25+l) 

^ Sr 



X [Esr - {E)feM~Esr/kBT) (13) 

with Cs(T) given by the sum over r for fixed S. The 
Cs{T) contributions in Eq. [13] are manifestly positive 
and are shown below. We also decomposed C(T) based 
on {E{T)). The results are similar at low T, the region 
of interest. 

Fig. [H] shows C(r, x) a.t X = 0.25 and 0.40 as a 
function of kBT/Ji for N — 18. The contributions of 
5 = 0, 1 and of 2 < 5 < 9 are indicated with dashed 
lines. Large E„i at a; = 0.40 gives a C{T) peak due 
to S = 0, 1 and a shoulder at higher T for S* > 2 
contributions that start with Eq w '2Em- Small Em 
at a; = 0.25 gives a single broad C(T) peak whose 




FIG. 12: Molar specific heat C{T) in Eq. dH in the BOW 
phase of H{x) at x = 0.25 and 0.40 for iV = 18 and PBC. 
The dotted and dashed lines are contributions in S = 0,1 and 
5 > 2, respectively, according to Eq. 1131 The inset shows the 
S = 0,1 contribution for x = 0.40 at iV = 16 and 18. 



Miitter and Wielath [36| reported C(T) results for 
8 < iV < 16 using a different numerical procedure and 
without resolving contributions in S. Their C'(T) and 
Xm{T) curves are quite similar to ours, but not identi- 
cal. In particular, their C(T, 1 /3) develops a shallow min- 
imum at = 16 that we do not see at either iV = 16 or 
18. Without proposing an explanation for the C{T, 1/3) 
maxima, Miitter and Wielath [s^ interpreted the dis- 
continuity of the C (T, x) maximum at the MG point as 
a transition from a "dimer" to a "frustrated" phase that, 
moreover, survived in the limit of large N. On the con- 
trary, our results indicate a single phase for a; > a;i and 
we understand the C(T) curves in Fig. [T^] in terms of 
5 = 0, 1 and S >2 contributions. We turn next to C(T) 
at large A'" = 2n. The fraction of singlets among the 2^" 
spin states is 



,/o(2n) = 



(2n)!2 



-2n 



(n+1)- 



n\{n + 1)! 



(14) 



The second expression follows from Stirling's approxima- 
tion and is accurate to within a few percent for 2n — 16 
or 18. The triplet fraction is 9n/o/(n -I- 2). Since /o de- 
creases as {2/N)-^/'^, the C(T) contribution from 5 = 
and 1 becomes negligible compared to 5 > 2 in the ther- 
modynamic limit. The inset in Fig. [T^] compares the 
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S' = 0, 1 maxima for N = 16 and 18. The N = 18 maxi- 
mum is lower as expected for Eq. [T31 The DMRG results 
of Feiguin and White [H for C{T, 1/3) at TV = 32 and 
64 show a smaU shoulder at low T that they attribute to 
finite-size effects. 



IV. DISCUSSION 

We have characterized the BOW phase of the linear 
spin-1/2 chain, H(x) in Eq. [U with frustrated Heisen- 
berg AF exchange Ji = J{1 — x) between neighbors 
and J2 = Jx between second neighbors. Exact VB 
methods yield the energies and eigenstates of finite 
systems with periodic or open boundary conditions. 
The BOW phase for a; > xi = 0.1943 has a broken 
Ci symmetry and finite magnetic gap Em{x). Our 
results are most accurate for 0.3 < x < 0.5 where large 
Em{x) ensures small finite-size corrections. Larger N, 
DMRG or other methods will be needed to character- 
ize the BOW phase with small Em{x) or B{x) at a; > 0.5. 

The BOW phase is particularly simple at a; = 1/3, 
the MG point, where the exact gs for PBC is either 
Kekule diagram \K1) or \K2) in Fig. [1] and \K1) for 
OBC. The BOW amplitude is B{x) in Eq. HI with 
-6(1/3) = 3/8 and a broad maximum in Fig. |6]at lower 
X than the Em{x) maximum in Fig. [5] Topological 
spin-1/2 solitons that reverse the bond order are the 
elementary excitations of the BOW phase, as shown in 
Fig. [Tlandilat x = 1/3 for odd N and OBC. The energy 
2Ew of two solitons corresponds for finite N to Em for 
parallel spins or to E^ for paired spins. We have also 
found the consequences of finite Em{x) on the molar 
spin susceptibility xm{T) and specific heat C{T). 

The magnetic properties of the EHM with parameters 
[/, V and t are closely related to H{x) when t >• [U — V). 
Van Dongen 37| mapped the EHM in the spin sector to 
H{x) with 



J(l -x) = Ji = 

xJ = J2 = 



4i^ 



U-V 
4t4 

{u - vY 



4J2 



(15) 



The HAF is the familiar limit t < (U-V) leading to 
X = 0. Increasing V at constant t and U amounts to 
increasing Jij J\ = xjiX — x). The CDW transition of 
the EHM is close to ^ = U 12. Since a continuous CDW 
transition 0, |1] requires t > U/7, the t <^ {U - V) ap- 
proximation fails at the BOW boundary of the EHM. The 
next term goes as and in addition to Ji , J2 contribu- 
tions, it adds [38| a four-spin contribution that requires 
going beyond H{x). Charge degrees of freedom cannot 
be neglected in the BOW phase of the EHM or of related 
models with Coulomb interactions. 



It is nevertheless attractive to approximate magnetic 
properties of Hubbard-type BOWs with H{x), much as 
HAFs have been used for Hubbard models. Finite Em 
implies an exponentially small xm{T) at low T followed 
by a roughly linear increase up to 0.15 in reduced units. 
As seen Figs. [10] and [TTl the xm maximum depends 
weakly on J2. Linear xm( T) vs T behavior is distinctly 
different from an HAF [3a| or a Hubbard model ^2^] with 
Em = and finite xm(0) that exceeds 60% of the xm 
maximum. Linear xm{T) following an onset has been ob- 
served [40] in several alkali-TCNQ sahs up to T w 450/^, 
the limit of their thermal stability. In our opinion, such 
Xm{T) in regular arrays are signatures of BOW phases 
with Em > in Hubbard-type models as well as in H{x). 

There are several reasons for considering 1:1 alkali- 
TCNQ salts as possible BOW systems. The strongest 
case is for Rb-TCNQ(II): its lOOK structure has 
regular stacks of TCNQ^ at inversion centers, negligible 
Xm{T) us to 150K and infrared spectra that indicate 
broken electronic inversion symmetry. Hubbard-type 
models have long been used for the magnetic, optical 
and electrical properties of quasi- ID organic ion-radical 
crystals [H [sj]. The singly occupied MOs of TCNQ" 
form a half-filled band. The BOW phase of Hubbard- 
type models is narrow, close to the CDW transition. 
1:1 alkali-TCNQ crystals are close [4^ to the CDW 
transition based on their electrostatic (Madelung) 
energy and the electronic structure of TCNQ~. The 
spin susceptibility of H{x) is encouraging for a BOW 
interpretation. More quantitative modeling will require 
values for t, U, V and other microscopic parameters. 

In summary, we have characterized the BOW phase 
of the linear spin-1/2 chain H{x) with frustrated first 
and second neighbor exchange. We exploited the exact 
gs at X = 1/3 for finite N to obtain the BOW ampli- 
tude 5(1/3) = 3/8, the magnetic gap Em, the spectrum 
of low-energy excitations, and topological spin solitons 
between BOWs with opposite phases. The spin chain 
makes possible a detailed examination of a BOW phase. 
While quantitative results are limited to H{x), the conse- 
quences of broken inversion symmetry, gs degeneracy and 
finite Em hold for BOW phases in general. The spin sus- 
ceptibility of H{x) is consistent with the unusual Xm{T) 
of alkali-TCNQ salts with regular stacks and provides ad- 
ditional support for the hypothesis [4l| that these salts 
are physical realizations of BOW phases. 
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